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We present an analytic study of a periodic Josephson array with long-range interactions in a 
transverse magnetic held. We find that this system exhibits a first-order transition into a phase 
characterized by an extensive number of states separated by barriers that scale with the system size; 
the associated discontinuity is small in the limit of weak applied field, thus permitting an explicit 
analysis in this regime. 



The essential features underlying the physics of glass 
formation are still unclear. In particular the possibility 
of glassiness in the absence of disorder is a question that 
has stimulated much recent activity pf. In this paper 
we present and study a periodic Josephson array with 
long-range interactions frustrated by a transverse mag- 
netic field (H). We show that it displays a first-order 
transition into a low-temperature phase characterized by 
an extensive number of states; these metastable solutions 
are separated by barriers that scale with the system size. 
The discontinuity at the transition is small in the limit of 
weak applied field, so that an explicit analytic treatment 
is possible in this regime. 
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Fig 1. A schematic of the array discussed in the text; 
the solid (shaded) lines are the horizontal (vertical) su- 
perconducting wires that are coupled by Josephson junc- 
tions to the vertical (horizontal) wires at each node. 



More specifically, the proposed array is a stack of 
two mutually perpendicular sets of N parallel wires 
with Josephson junctions at each node (Figure 1) that 
is placed in an external tranverse field. The classical 
thermodynamic variables of this system are the super- 
conducting phases associated with each wire. Here we 
shall assume that the Josephson couplings are sufficiently 
small so that the induced fields are negligible in compar- 
ison with H. We can therefore describe the array by the 
Hamiltonian 
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where J m n is the coupling matrix 



J = 



J 
Jt 



(2) 



with J jk = -j2^exp(2niajk/N) and Sj(s m ) is a N(2N)- 
element vector whose components are exp(i0 m ), where 
4> m are the superconducting phases of the wires. Here we 
have introduced the flux per unit strip, a — N HI 2 / <po, 
where I is the inter-node spacing, 4>q is the flux quantum, 
and the normalization has been chosen so that T c ~ Jo- 
Because every horizontal (vertical) wire is linked to 
every vertical (horizontal) wire, the number of nearest 
neighbors in this model z = N; we can therefore study 
it with a mean-field approach. For the same reason the 
free energy barriers separating its low T solutions scale 
with N, in marked contrast to the situation in a conven- 
tional 2D array H; thus each of these solutions results 
in a metastable state. If the field is such that the num- 
ber of solutions is extensive (i.e. the phase-ordering unit 
cell scales with the system size) this array displays glassy 
behavior in the absence of disorder. 
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A similar long-range network with areal disorder was 
studied previously |; for a > 1/iV the system under- 
goes a spin glass transition which was mapped onto the 
Sherrington-Kirkpatrick model for the case a»l for all 
T. Physically this behavior occurs because the phase 
differences associated with the couplings, Jjk, acquire 
random values and fill the interval (0, 2ir) uniformly for 
a ^> 1/N. In the absence of randomness these phases still 
cover the interval for a ^S> 1/N, so that one expects that 
glassy behavior survive Q . Naturally there will be spe- 
cial values of the field for which the number of solutions 
is finite. However there is no such commensurability for 
1/N <C ol < 1; e.g. the unit cell corresponding to a = 1 
consists of all the wires traversing a strip and thus has 
sides of order 0(N)l. 

The cumulants of the coupling constants display a 
more quantitative similarity between the random and the 
periodic arrays. In both cases 
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which can be understood physically since only horozontal 
(vertical) wires separated by d < l/a contribute coher- 
ently to this sum independent of disorder. 

We now begin a more quantitative analysis of this pe- 
riodic array using a modified Thoulcss- Anderson-Palmer 
approach || . Because the interactions are long-range, we 
can integrate out the thermal fluctuations and define a 
"site magnetization" m, = {z{)t- The free energy is 



F{m k } = - m%Jkjrnj -T^S(m fc ) + F {m k } 



-^2{h k m* k + h* k m k ) 



S{ m ) = S -\m\ 2 ~-\m\ 4 + O(\m\ 6 ) 



(4) 



where J k j is defined in 



Here we have introduced 



hi, a field conjugate to and the single site entropy 
S{m); F {m k } is the Onsager reaction term. At high 
temperatures F {m k } ~ 1/T and thus is small; fur- 
thermore, as we show below, the transition occurs at 
T c « Jo/ \fa (as in the disordered case) so that we ex- 
pect the feedback effects to be perturbative for small a 
(1/N <C a <C 1). Thus we shall ignore F {m k ] in our 
initial study of F{m k }, but will consider it later. 

In the high temperature phase m, = 0. This "param- 
agnetic" state becomes unstable when the quadratic part 
of F{m k } acquires negative eigenvalues, i.e. at T < X max 
where \ max is the largest eigenvalue of J. Due to the 
structure of J (see @), any eigenvalue A of J.P results 
in two eigenvalues A = ±\/A of J, a property we exploit 
for convenience. In the limit N — > oo 
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which can be diagonalized by a Fourier transformation 
to yield 
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Therefore the largest eigenvalue of J is X ma x — Jo/y/a; 
the corresponding eigenfunctions associated with one set 
of parallel wires are plane waves with momenta in the 
interval < p < 2ira so that the degeneracy of this 
eigenvalue is aN 0. 

Therefore, in the absence of the feedback effects, the 
transition occurs at T c q = Jo/y/a into a phase that is 
characterized by a linear combination of the aN states 
described above. Thus the number of possible metastable 
configurations just below the transition should be exten- 
sive, indicating the onset of glassiness 0. This should 
be contrasted with the situation in the Sherrington- 
Kirkpatrick model where, in the thermodynamic limit, 
the highest eigenvalue of the exchange matrix is non- 
degenerate and the metastability measure of the or- 
der parameter (plateau-onset xq in Parisi function q(x)) 
grows continuously beneath the transition. More gener- 
ally, in disordered spin glasses a sudden appearance of 
extensive metastability implies a first-order jump in the 
order parameter B. Unfortunately in the periodic case 
we have not yet determined the analogous order param- 
eter; however we show below that the Onsager feedback 
term drives the transition weakly first-order, consistent 
with our expectations. 

We use the locator expansion to determine the lead- 
ing order contributions in rm to F {m k }, it is based on 
the expression for susceptibility 
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A- J 



where A is a diagonal "locator" matrix. This form of 
susceptibility is a consequence of large z; to prove it one 
can use the high temperature expansion and show that 
all diagrams renormalizing the J term are small in 1/z. 
Once we have solved for A, we can reconstruct the free 
energy from 
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that gives us the free energy which contains both the 
entropy term ^ TS(mj) and the Onsager term F {rrii}. 
Although it is possible, in principle, to determine the el- 
ements of A ab initio it is easier to do so using the iden- 
tity Xjj = 1 ~^" J ^ which follows from the fluctuation- 
dissipation theorem. We expand the expression for the 
susceptibility (^) and, due to the structure of J , find 
that only terms even in J contribute. Collecting and 
resumming these terms, we obtain 
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We first consider the leading-order term (setting m,j = 0) 
which will yield the quadratic part of F{rrij}. To this or- 
der A is independent of the site label {A — AY), so it is 
convenient to solve this equation by Fourier transforma- 
tion and, using (JsJ) , we obtain 
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The zeroth-order term in a yields A = T, consistent 
with our assumption above that the reaction term is small 
in a. 
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Fig 2. Plot of the locator A(T) obtained as a numerical 
solution of Eq. (10) for Jq = I, a = 0.1; the horizontal 



line denotes the critical value A(T) = X r , 
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The solution A(T) has two branches A ± (T) displayed 
in Figure 2. Since A(T) is always larger than Jo/^/a 
on the upper branch (A + (T)), within this approxima- 
tion the high temperature phase is locally stable down to 
zero temperature. It is therefore impossible to reach the 
low temperature phase from high temperatures without a 
discontinuity in A(T). A lower limit on the discontinuity 
is the gap between two branches 5 A — | y/aJo <C Jo / \fa 
(recall that we are taking a -C 1). We note that although 
the feedback terms are generally small, they become im- 
portant near the transition and here change its order. We 
can estimate the corresponding minimal jump in \mi\ 2 



neglecting the feedback effects on the quartic terms (i.e. 



using = [A(T) - A r , 



= iTK| 4 ), to be 
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= a. In the vicinity of T c0 = X max the disconti- 



nuity becomes larger (AA = X max - A (T c0 ) = J /V2); 
since we expect the transition to occur in the vicinity 
of T c o, the most likely value of magnetization jump is 
A | mj | 2 = \^2a, as is also borne out by a more sophisti- 
cated calculation Jl0[ . 

In the regime of strong frustration, when a ~ 1, the 
transition becomes strongly first order. The case a = 1 
has a special interest because in this situation J becomes 
a unitary matrix and cumulants of J resemble those of 
the SK model: ^TrJ" 2 " = J o"- As a result, the calcula- 
tion becomes simpler. The leading order locator equation 
( |l0| ) can be solved explicitly 
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suggestive of a strong first order transition around T ~ 
Jo- Certainly, the existence of such a transition cannot 
be determined by a m-expansion. Instead, we adopt a 
variational approach and demonstrate that there exists 
a "magnetic" phase whose free energy is lower than the 
"paramagnetic" free energy F + — —2NTSo(T) associ- 
ated with the upper branch A + (T). We choose the trial 
state 
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exp(-«— 7J ) 



TO ^ = W V2 exp(j - 7 -i fc2) 

where 7 is a variational parameter; since \m,j 
free energy of this state is equal to its energy 
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where the dimensionless function £(7) satisfies £(7) = 
£(1/7) and £(7) = 7 at 7 < 1. Clearly at T < J this 
state has a lower free energy than the paramagnetic state 
associated with A + (T), and thus we conclude that it is 
unstable at low temperatures. 

We emphasize that it is difficult to observe a param- 
agnetic instability in the free energy expanded in m.;; we 
have already seen that its quadratic part remains posi- 
tive definite on the upper branch A + (T). Now we show 
that this conclusion does not change when quartic terms 
are taken into account. To do so we assume the most 
general form of the free energy possible at a = 1 

F {mj} = Y^mUhju, + Y}A{T)\m 3 \ 2 + ~&Kf] 

kj j 

-^2(h*mj + c.c.) + ^^2a k3 \m k \ 2 \m 3 \ 2 (15) 
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where the coefficients dkj 



if {k,j} label parallel 



wires and a k j = a- otherwise. We solve the equation 
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ohtOhjdh'. ' [t also satisfies Xkj 
Comparing these two expressions we determine the coef- 
ficients a + , a_ and b: 



b = T 
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The quartic term of the free energy is minimized when 
|mj| 2 = m 2 , but still remains positive for the upper 
branch (A+(T) > J ) 



F(m) = 2N 



(A(T) - 4W + \T§^ 



T 2 
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Of course, a glass transition is best characterized by 
its dynamical features, e.g. relaxation spectrum, hys- 
teresis, ageing and memory effects; we have not yet in- 
vestigated these properties of the periodic array discussed 
here. Hopefully the order parameter associated with this 
glass transition will also emerge from such a study, since 
we expect it to be a dynamical quantity. The thermody- 
namic behavior of a similar array with a reduced number 
of neighbors is another open question which could be 
important in understanding the relevance of long-range 
models to experimental glasses. 
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implying that all solutions remain stable to quartic order. 

In summary, we have presented a non-random model 
which exhibits a first-order transition to a phase where 
the number of metastable states is extensive, consistent 
with one's general expectations for glassiness in the ab- 
sence of disorder 0] . In this array the interaction between 
wires is long-range, and thus the barriers separating the 
low-temperature metastable states scale with the size of 
the system. Since the number of nearest-neighbors z is 
large we could perform an initial investigation using a 
modified TAP approach; furthermore the glass transition 
is weakly first-order in the limit of small a (| « a « 1), 
the flux per unit strip, so that it could be studied ana- 
lytically. 

The results presented here can be tested in numerical 
simulations. The system will have a signatory superfluid 
stiffness, p s , which can be measured by applying a twist 
$ to the the horizontal(vertical) wires at their respective 
boundaries |i~l| ; we predict that p s = will appear dis- 
continuously at the transition T c (a) and that the jump 
will scale with a. Thermal cycling should indicate an ex- 
tensive number of metastable states, thereby confirming 
the low-temperature glassy phase predicted here. 

This transition could also be probed in fabricated ar- 
rays with negligible induced fields; more specifically we 
require 
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N 5 / 2 l 
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which constrains the size of the array |12|. Here T c (a) 
should be observable by resistive measurements. Alter- 
natively one can infer the superfluid stiffness from induc- 
tive experiments. The critical current below T c (a) should 
jump as well. 
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